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Abstract 



> 

We calculate the 1 Sq pairing gap in nuclear matter by adopting the "in- medium 
Bonn potential" proposed by Rapp et al. [e-print nucl-th/9706006], which takes into 

OO '• 

. account the in-medium meson mass decrease, as the particle-particle interaction in 



the gap equation. The resulting gap is significantly reduced in comparison with the 
one obtained by adopting the original Bonn potential. 
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Superfluidity in infinite hadronic matter has long been studied mainly in neu- 
tron matter from a viewpoint of neutron-star physics such as its cooling rates. 
As a way of description, relativistic models are attracting attention in addi- 
tion to traditional non-relativistic nuclear many-body theories. Since Chin and 
Walecka succeeded in reproducing the saturation property of symmetric nu- 
clear matter within the mean-field theory (MFT) [1], quantum hadrodynamics 
(QHD) which is an effective field theory of hadronic degrees of freedom has 
described not only infinite matter but also finite spherical, deformed and ro- 
tating nuclei successfully with various approximations [2,3]. These successes 
indicate that the particle-hole (p-h) channel in QHD is realistic. In contrast, 
relativistic nuclear structure calculations with pairing done so far have been 
using particle-particle (p-p) interactions borrowed from non-relativistic models 
such as Gogny force. Aside from practical successes of this kind of calculations, 
the p-p channel in QHD itself is a big subject which has just been started to 
study. 

The first study of this direction was done by Kucharek and Ring [4]. They 
adopted, as the particle-particle interaction (f pp ) in the gap equation, a one- 
boson-exchange interaction with the ordinary relativistic MFT parameters, 
which gave the saturation under the no-sea approximation. The resulting max- 
imum gap was about three times larger than the accepted values in the non- 
relativistic calculations. [5-9]. Various modifications to improve this result 
were proposed [10-12] but this has still been an open problem. One of such 
modifications is to include the p-h and the N-N polarizations, the former of 
which is effective for reducing the gap in the non-relativistic models [13,14]. 

From a different viewpoint, Rummel and Ring [15,2] adopted the Bonn poten- 
tial [16], which was constructed so as to reproduce the phase shifts of nucleon- 
nucleon scattering in free space, as v pp and obtained pairing gaps consistent 
with the non-relativistic studies. Since the single particle states are determined 
by the MFT also in this calculation, explicit consistency between the p-h and 
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the p-p channels is abandoned. However, assuming that the MFT simulates the 
Dirac-Brueckner-Hartree-Fock (DBHF) calculation using the Bonn potential, 
we can consider that the consistency still holds implicitly. A possible extension 
of this study is to take into account the in-medium changes of the properties of 
the mesons which mediate the inter-nucleon forces. Among them, the change 
of the mass, that is, the momentum-independent part of the self-energy, has 
been discussed in terms of the Brown-Rho (BR) scaling [17], the QCD sum 
rules [18], and some other models both in the quark level and in the hadron 
level [19,20]. Although there still has been theoretical controversy [21], some 
experiments look to support the vector meson mass decrease [22]. In addition, 
the momentum-dependent part also attracts attention recently [23-26] . In the 
present work, we assume only the momentum-independent vector meson mass 
decrease. A simple way to take into account this meson mass decrease in the 
meson exchange interactions is to assume the BR scaling. Actually Rapp et 
al. [27] showed, based on this, that the saturation property of symmetric nu- 
clear matter and the mass decrease of the vector mesons were compatible. 
So we adopt their "in-medium Bonn potential" as v pp in the gap equation. 
Since the gap equation takes the form such that the short range correlation 
is involved [28,29,7], we use this interaction in the gap equation without the 
iteration to construct the G-matrix. 

As described in ref. [4], meson fields also have to be treated dynamically 
beyond the MFT to incorporate the pairing field via the anomalous (Gor'kov) 
Green's functions [30]. The resulting Dirac-Hartree-Fock-Bogoliubov equation 
reduces to the ordinary BCS equation in the infinite matter case. We start 
from a model Lagrangian for the nucleon, the a boson, and the uj meson, 

C = $(i lll d>* — M)ip 
+ \(dr*)(d»v)-\mlo* 
+ g^aij) - g^^ip, 
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The actual task is to solve the coupled equations for the effective nucleon mass 
and the 1 S pairing gap: 



where v pp (p,k) is an anti-symmetrized matrix element of the adopted p-p 
interaction, 



with an instantaneous approximation (energy transfer = 0) and an integration 
with respect to the angle between p and k to project out the S'-wave com- 
ponent, and (u>°) is determined by the baryon density. Here we consider only 
symmetric nuclear matter (7 = 4) and neutron matter (7 = 2). Conceptually 
the numerical integrations run to infinity in the present model where the finite- 
ness of the hadron size is considered in v pp . Numerically, however, we introduce 
a cut-off A c ; its value is chosen to be 20 fm -1 so that the integrations converge. 
If we adopt the Bonn potential as v pp and include the p meson and the cubic 
and quartic self-interactions of a with choosing the NL1 parameter set [31] for 
the mean field, these equations reproduce the results of ref. [15]. Here we note 
that there are pros and cons as for the self- interaction terms [32,33]. First we 
developed a computer code to calculate the Bonn potential for this channel 
from scratch, and later confirmed that our code reproduced outputs of the one 




e fc = v / k 2 + M* 2 + ^(^°), 



(2) 



v PP (p,k) 

= (ps',ps'|i;pp|ks,ks) - (ps', ps'\v pp \ks, ks), 



(3) 



in ref. [34]. 
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Then in the present work we adopt the in-medium Bonn potential of Rapp et 
al. [27] although this is, as yet, unpublished. To construct this potential, they 
applied the BR scaling after replacing the a boson in the original Bonn-B po- 
tential with the correlated and the uncorrelated 2n exchange processes. Then 
they parameterized the obtained nucleon-nucleon interaction by three scalar 
bosons 0i, o" 2 and "rest-a", in addition to n, r], p, u and 5. Here <7\ and a 2 with 
density-dependent masses and coupling constants simulate the correlated 2n 
processes and the "rest-a" simulates the uncorrelated ones. Therefore the ac- 
tual tasks are adding two extra scalar bosons to the original Bonn-B potential 
and applying the BR scaling, 



where po is the normal nuclear matter density, A PiW are the cut-off masses in 
the form factors of the meson-nucleon vertices, which is related to the hadron 
size. The scaling parameter C is chosen to be 0.15. Aside from a slight tuning 
of the coupling constant of 5, other parameters are the same as those in the 
original Bonn-B potential. All the potential parameters are presented in Tables 
I and II in ref. [27]. Note that this scaling applies only to the quantities in i> pp , 
since the MFT with its effective nucleon mass guarantees the saturation and 
we confirmed that the effects of pairing on the bulk properties are negligible 
except at very low densities in the present model. As for the 7r-N coupling in 
the potential, we adopt the pseudovector one conforming to the suggestion of 
chiral symmetry [35] and actually the pseudoscalar one in the Bonn potential 
was shown to lead to unrealistically attractive contributions in the DBHF 
calculation [36]. 

The result is presented in Fig.l. This shows that applying the BR scaling re- 
duces the gap significantly in comparison with the original Bonn-B potential. 
The potentials themselves are given in Fig. 2. The relation between the reduc- 
tion of the gap and the change in the potential, i.e., the shift to lower momenta, 



M* 



m* 




A* 



M 



(4) 



m. 
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can be understood as follows: As discussed in refs. [6,15], A(k) determined by 
the second equation of (2) exhibits a nodal structure similar to v pp (k F ,k) as 
functions of k with the opposite sign and some possible deviation of zeros. Con- 
sequently both the low- momentum attractive part where A(k) > and the 
high-momentum repulsive part where A(k) < give positive contributions to 
A(fcp), and Fig. 2 shows that the both of these contributions are reduced in the 
case of the in-medium Bonn potential. This is the main reason why A(k F ) is 
reduced by applying the BR scaling. Actually we confirmed that this is mainly 
due to the mass decrease of the vector mesons, not of the nucleon, as shown in 
Fig. 3 although also the latter produces some reduction of the gap. Note here 
that the saturation in the DBHF calculation is guaranteed only when both M 
and scaled. The in-medium meson mass decrease is brought about 

by the N-N polarization in hadronic models. Although its explicit inclusion in 
the gap equation has not been reported yet, it is conjectured in ref. [12] that 
its inclusion simultaneously with the p-h one will reduce the gap as that of 
the p-h one in the non-relativistic models [13,14]. In this sense, the present 
result that applying the BR scaling reduces the gap is consistent with this 
conjecture, although the reduction in refs. [13,14] is more drastic than that in 
the present work. An additional element is that the cut-off masses in the form 
factors of the meson-nucleon vertices are also scaled. Since they appear in the 
form 

where q is the momentum transfer [16], simultaneous reductions of A a and 
m a (a = p,co), lead to a further reduction of the repulsion due to the vector 
mesons at large- |q|. Since the higher the background density becomes the more 
nucleons feel the high- momentum part of v pp , the in-medium reduction of A p w 
leads to an additional reduction of the gap in the high-density region. 

Finally, use of the NL1 parameter set for the mean field brings negligible 
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changes and the result for neutron matter is very similar to those presented 
here for symmetric nuclear matter. 

To summarize, we adopted the in-medium Bonn potential, proposed by Rapp 
et al., which takes into account the in- medium meson mass decrease in a sim- 
ple form and is compatible with the nuclear matter saturation in the DBHF 
calculation, as the particle-particle interaction in the gap equation. The result- 
ing pairing gap in nuclear matter is significantly reduced in comparison with 
the one obtained by adopting the original Bonn potential. Here we should 
note that the uncertainty in the gap value deduced from various studies is 
still larger than the strength of the reduction found in the present work, and 
therefore, further investigations are surely necessary both relativistically and 
non-relativistically. 
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Figure Captions 



Fig.l. Pairing gap in symmetric nuclear matter at the Fermi surface as func- 
tions of the Fermi momentum. Dotted and solid lines indicate the results 
obtained by adopting the Bonn-B and the in-medium Bonn potentials, re- 
spectively. Note that accuracy is somewhat less around A(fcp) ~ in our 
code. 

Fig. 2. Matrix element v pv (k-p, k) as functions of the momentum k, with a Fermi 
momentum k F = 0.9 fm _1 . Dotted and solid lines indicate the results obtained 
by adopting the Bonn-B and the in-medium Bonn potentials, respectively. 

Fig. 3. Pairing gap in symmetric nuclear matter at the Fermi surface as func- 
tions of the Fermi momentum. Solid, dotted and dashed lines indicate that 
the results obtained by reducing according to Eq. (4) only the nucleon mass, 
only the vector meson masses, and the both, respectively. 
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Pairing Gap at the Fermi Surface 

parameter set a-co, two types of Bonn are used in the gap equation 
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Pairing Gap at the Fermi Surface 

parameter set a-co, In-Medium Bonn is used in the gap equation 
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